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A linearly coupled chain of spin-polarized quantum dots is investigated under the condition that 
the number of electrons is equal to or less than the number of the dots. The chemical potential of 
the system, ^jv = E{N) - E{N - 1), satisfies, (/ijv + ^i]Vf+2-jv)/2 ~V-\-2t {N, Ne, V, E{N) and 
t are the number of electrons, the number of dots, and the strength of nearest neighbor electron- 
electron interactions, the total groundstate energy and the hopping integral between two adjacent 
dots). This property will be reflected in the spacing between the conductance peaks. The electron 
density structures are determined using a quantum Monte Carlo method. As the number of electrons 
is varied several correlated structures are found that are commensurate/incommensurate with the 
periodic dot system. Hartree-Fock theory fails to predict the correct electronic structures of this 
system because several nearly degenerate solutions exist. 



PACS numbers; 73.20.Dx, 73.20.Mf 



Recent advances in nano-fabrication techniques have 
made it possible to make quantum dots 0J^. These 
structures have several similarities to atoms and are 
called artificial atoms. The physical parameters of these 
systems may be controlled easily: the electron density 
can be varied significantly by the substrate voltage Vg 
and the range of electron-electron interactions can be 
controlled by changing the distance to the metallic layer. 
When these artificial atoms are coupled structures similar 
to a molecule may arise. Recently numerous groups have 
started to investigate whether coupled two-dot systems 
really have molecular properties. Properties of a two-dot 
system have been studied both experimentally and 
theoretically |^,^ . When many dots are coupled an artifi- 
cial one-dimensional crystal may be created. It is possible 
to couple up to ten or more dots together since the size 
of the total system can be made smaller than the phase 
coherence length. Such a system has energy bands sep- 
arated by energy gaps, and transport properties of such 
a periodic crystal have been investigated experimentally 

When the length of such a periodic crystal is infi- 
nite the system exhibts many intriguing properties. De- 
pending of the values of the physical parameters the 
system can have properties of a Luttinger liquid or a 
generalized Wigner crystal It is unclear whether 

these properties remain in a finite system of coupled 
dots. The purpose of this paper is to examine general 
properties of a linear chain of spin-polarized dots un- 
der the condition N < N^. We find the following prop- 
erties. The chemical potential of the system satisfies, 
{fiN + MWf+2-Af)/2 « V -\-2t. This property of the chem- 
ical potential is reflected in the spacing between the con- 



ductance peaks Q . The electron density structures, de- 
termined by a quantum Monte Carlo method, show sev- 
eral types of commensurate/incommensurate structures 
as the number of electrons is varied. The groundstate 
electronic properties of these structures cannot be de- 
scribed by mean field theory (Hartree-Fock theory) ex- 
cept for a few isolated cases. The reason for this failure 
is the existence of many nearly degenerate solutions. 

As shown by recent experiments when N < Ni the sys- 
tem can be easily spin-polarized [10|. The Hamiltonian 
of such a coupled dots is given by |13| 

H = Ht + Hv 

= ('^tic. + c+Q+i) + 2tN (1) 

1=1 

Ne-l 



The operator Ci(cf) destroys (creates) an electron in the 
i-th dot. The term 2tN allows to measure single particle 
energies form zero. 

The presence of the end sites at 1 and Ni breaks 
particle-hole symmetry. However, the system can be 
mapped into a ring of (A^^ + l)-coupled dots by imposing 
the periodic boundary condition CArj+2 = ci and requir- 
ing the electron occupation number at the site A'^ -|- 1 to 
be zero: 



Ne + l 
i=l 



Q + 



etc,. 



+V ^ niUi+i + WNi+iriNi+i + "itN. (2) 

i=l 

When the strength of the impurity potential energy, 
Wlq^+l^ is large and positive the occupation number at 
the site iV^ + 1 will be zero. We rewrite the Hamilto- 
nian using the hole creation and hole occupation number 
operators = q and rrii — 1 — rii: 

H = H' + {V + 2t){2N - Ni^l), (3) 

where 

Ne + l 

i=l 
Ni + 1 

+V ^ minii+i 
1=1 

+WN,+i{l - mN,+i) + 2tNh. (4) 

The quantity Nh is the total number of holes, equal to 
Ni + 1 — N . From this transformation we find the relation 

E{N) = E'irn + l-N) + {V + 2t){2N - Ne - 1), (5) 

where E{N) is the groundstate energy of H with N elec- 
trons, and E'{N) denotes the counterpart of H' . For 
sufficiently large Nf the impurity contributions to the 
groundstate energies are negligible in both H and H' , 
and we may set 

E'{N) « E{N) (6) 
^From this it follows 

{^lN + ^^N,+2-N)/2^v + 2t. (7) 

For comparison it should be noted that complete 
electron-hole symmetry of a ring of Ne dots leads to 

{^lN + ^^N,+l-N)/2^v + 2t. (8) 

Note that the second subscript in Eq. (^) contains f in- 
stead of 2. In a finite linear chain or in the presence of 
disorder this relation is not valid because electron-hole 
symmetry is broken. 

We have solved Eq. (I) by a Hartree-Fock method. The 
ground state is determined by the following equations 

HHF\k) = ek\k) (9) 

Here the Hartree-Fock Hamiltonian is given by 

Hhf = -tYl (clc,+i + + ^Xiclc, 

i=l 1=1 
Ne-1 

(10) 

1=1 



where 

X, = y((n,+i> + 

The electron density {{rii)} and the Hartree-Fock eigen- 
statcs {|fc)} should be determined self-consistently. Note 
that (no) — {riMi+i) = 0. The total energy of iV-electron 
system is given in terms of the Hartree-Fock eigenstates 
{\k)}: 

N N 

E{N) = E^fe - 2 ^ {{kq\Hv\kq) - {kq\Hv\qk)) , 

k=l k,q=l 

(11) 

where 

ATf-l 

{kq\Hv\kq) ^ 2V ^ \a,{k)\''\a,+i{q)\'' , 
1=1 

Ne-l 

{kq\Hv\qk) ^ 2V ^ a^{k)* a^{q)a^+i{k)* ai+i{q). 

The amplitude of the fcth eigenstate at site i is ai{k). 

Fig. 1 displays the lowest HF solutions for an odd value 
of Ni = 7. Since they have nearly degenerate ener- 
gies we must include quantum fluctuations to find the 
true groundstate. Fig. 2 displays the density profiles ob- 
tained using a quantum Monte Carlo method jlj]. We 
used typically 500 000 - 1 000 000 Monte Carlo steps 
to measure the values and the statistical error is about 
the size of the symobols. Quantum Monte Carlo and HF 
results differ significantly. The HF theory overestimates 
the repulsive interactions, and consequently favors struc- 
tures with more oscillations. A commensurate structure 
is found at = 4. Generally, when Ne is odd a com- 
mensurate state exists for N = [Ni -I- 1)/2. This type of 
states have interesting optical properties [p] . Quantum 
Monte Carlo results show that the number of peaks is 
equal to N when N < Ng/2. The number of minima is 
equal to A^ - TV for N > Ne/2. 

Fig. 3 displays quantum Monte Carlo density profiles 
for an even value of A^^ = 14. In both cases the following 
structures are present: (•, •, o, •, •, o), (o, o, •, o, •, o, o, •) 
where the symbols • and o denote relatively large and 
small occupation numbers. HF or classical theory fails 
to predict these structures. Note that when Ni/N > 1/2 
the occupation number increases as the site index moves 
away from the center and the values of Ui are large at 
1 = 1 and Ni. This is because the electrons feel strong 
mutual repulsion. The opposite is true for Ng/N < 1/2. 
Note that, in contrast to infinite systems, the ground- 
states do not resemble a liquid state even for V = 1. 
Generally when A^^ is even a commensurate (periodic) 
state is absent, since two classically degenerate state ex- 
ists. Again we note that the number of peaks is A^ when 



2 



N < Ni/2, and the number of minima is Ni ~ N when 
N > Ni/2. 

The inset in Fig. 4 displays N versus jiN for A'^ = 15 
and V = A. The distance between the {N + l)-th and 
iV-th peaks is equal to hn+i — fJ'N The constant 

charging model would fail to account for these results 
since the separations between the peaks are not a con- 
stant. We have tested the accuracy of Eq. (7) numerically 
by plotting {iJ,]\[+fi]y^+2-N)/2 as a function of N. We see 
in Fig. 5 that even for small Ni Eq. (7) is well satisfied. 

Kouwenhoven et al. pO| ] have shown that their exper- 
imental data are better described by assuming presence 
of some disorder. We have also carried out a similar cal- 
culation in the presence of an impurity, and find that our 
result is in qualitative agreement with that of Kouwen- 
hoven et al. It would be interesting to test the valid- 
ity of Eq. (5) experimentally in a clean coupled-dot sys- 
tem without impurities. Strong correlation effects should 
be more visible in conductance positions the larger the 
ratio V/t is. A finite chain of linearly coupled quan- 
tums dots is well suited for observing both commensu- 
rate/incommensurate electronic structures. 
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FIG. 1. Occupation numbers, n^, as a function of site index, i. The three lowest energy states of a Hartree-Fock calculation 
are plotted. The strength of nearest-neighbor electron-electron interactions, the electron number, and the number of dots are 
V = 5, N = 3, and Ne = 7. Here the energies are measured in units of t. The graphs are shifted vertically for the sake of 
clarity of display. 

FIG. 2. Quantum Monte Carlo results of occupation numbers for six different electron numbers A'^. The parameters are 
1/ = 5 and Af = 7. 



FIG. 3. Quantum Monte Carlo results of occupation numbers for nine different electron numbers N. The solid lines are for 
V — 4 while dotted lines are for V = 1. The total site number is Ne — 14. 
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FIG. 4. Chemical potential, hn, measured in unit oit is plotted as a function of A'^. Results for {Ni, V) equal to (15,4), (14,4), 
and (14,1) are represented by circles, inverted triangles, and triangles. The inset displays N vesus /itjv for (Ne,V) = (15,4). 
The energy difference between two adjacent values of /liw corresponds to the spacing between two nearby conductance peaks. 

FIG. 5. The quantity (jUAr+/xjv^+2-jv)/2, measured in unit of t is plotted as a function of N to test whether it is approximately 
equal to a constant V+2t. Results for {Ne, V) equal to (15,4), (14,4), (14,1), (7,4), and (7,1) are represented by circles, inverted 
triangles, and triangles, diamonds, and squares. 
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